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ABSTRACT 


Estimates of low flow and flood frequency in several southwestern Wisconsin 
basins have been significantly improved by determining land cover from Landsat 
imagery. With the use of new estimates of land cover in multiple-regression 
techniques, the standard error of estimate (SE) for the least annual 7-dav low 
flow for 2- and lO-yesir recurrence intervals of ungaged sites were lowered by 
9 percent each. The SE of flood frequency in the "Driftless Area" of Wisconsin 
for 1C-, 50- , and lOC-year recurrence intervals were lowered by ll< percent. 

Four of nine basin characteristics determined from satellite imagery were significant 
variables in the multiple-regression techniques, whereas only 1 of the 12 
characteristics determined from topographic maps was significant. The percentages 
of land-cover categories in each basin were determined by merging basin boundaries, 
digitized from quadrangles, with a classified Landsat scene* Both the basin 
boundary X-Y polygon coordinates and the satellite coordinates were converted to 
latitude-longitude for merging compatibility. 
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IMPROVING STREAt>IFLOW ESTIMATES THROUGH THE USE OF LANDSAT 


By G. J. Allord and F. L. Scarpace 
INTRODUCTION 

Comprehensive water-resojrces management needs information on streamflcv 
characteristics. Complete streamflow information requires techniques to extrapolate 
available information to ungaged sites. 

Thomas and Benson (1970) demonstrated that streamflow characteristics can 
be related to basin characteristics. This method regresses a streamflow characteristic 
against the basin characteristics of the gaged sites. The technique provides 
equations that define streamflow characteristics in terms of statistically 
signi'f’icant basin characteristics. These equations then may be used to estimate 
streamflow characteristics at ungaged sites where long-term records are not 
available. 

Land use is one of the basin characteristics shown to influence streamflow. 
Traditionally, the U.S. Geological Survey (USGS) has measured these land-use 
areas in Wisconsin from topographic maps, these being the best data available on 
a statewide basis. Landsat imagery can be used to provide up-to-date land-cover 
information over large areas with the advantage of special and temporal consistency. 

Two basins in southwestern Wisconsin were selected to demonstrate that 
determination of land use through Landsat imagery could improve the equations 
developed to predict streamflow characteristics at ungaged sites. The basins 
were selected because relations had been developed for the two basins, and for 
many subbasins within each basin, based on basin characteristics detennined by 
conventional techniques. 
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BASIN DESCRimON 


The tvo basins selected for study are the lower Wisconsin River basin and 
the Pecatonica-Sucar River basin (fig. l). The major part of both basins lies 
in the "Driftless Area" of Wisconsin. This area is termed "Driftless Area" 
because it is not covered by glacial drift, except for some alluvial and lake 
deposits in valleys (fig. l). The topography of the "Driftless Area" consists 
of rugged, steep-walled valleys and high relief. Agriculture, especially in the 
lower Wisconsin River basin, is confined to valley floors and hilltops because 
of the steepness of the slopes. The remaining part of each basin consists of 
rolling hills and wide valley floors fonned by the unconsolidated glacial deposits 
that overlie the sandstone bedrock. 

LANDSAT IMAGE CLASSIFICATION 

The Landsat scene (5386-15325, May 9, 1976) that covers the two river 
basins was analyzed to provide land -cover information. This particular Landsat 
scene was selected because it covered both basins, the season was well suited 
for determination of land-cover types, and atmospheric conditions and image 
quality were acceptable. 

The digital-image processing of the Landsat scene (fig. 2) was begun on 
3/ 

the Image-lOO System— at the EROS Data Center. Final classification, digitizing 
of the cartographic data base, and the merging of the map information with the 

3/ 

Landsat scene was completed on the University of Wisconsin-Madison Univac 1110=^ 
Computer System. The land-use and land-cover classification system suggested by 
Anderson and others (1976) was used as a guide for the definition of Level I 
(Landsat type of data) land -use categories. 


" The use of trade names in this report is for identification only and 
does not imply endorsement by the U.S. Geological Sui'vey. 
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The following; land-use categories were classified: forest, grassland, 

water, wetland, mixed vegetation types, bare soil, and cropland. Areas htiown to 
be representative of each land-use category were selected, and the spec ... i 
signatures from each area were determined. These s''.gnatures then were used to 
classify the rest of the image into land-use categories. Pinal classification 
of the land uses was done through the use of an unsupervised parallelepiped 
classifier with a mininuja distance-to-mean check for potentially overlapping 
classes. The results of this classification were stored on a computer tape and 
merged with the digitized basin boundaries provided by the technique described 
in the following section. 

CARTOGRAPHIC DATA BASE 

Basin boundaries, drawn on U.S. Geological Survey 1:2U,000 and 1:62,500 
quadrangles, were digitized at the University of Wisconsin-Madison Cartographic 
Laboratory. The line segments were converted from digitiser coordinates to 
latitude and longitude coordinates and stored in a data base, identified by the 
site-identification code assigned by the USGS downstream site-numbering system. 

The site-identification code is used to retrieve all the coordinates that would 
draw the polygon outline of the basin. The computer program provides the 
outline of a basin composed of several subbasins by compiling only the outside 
boundaries of the individual subbasins. 

MERGING LAI.’DSAT AND CARTOGRAPHIC INF0R?.IATI0N 
Although the geometry of the Landsat scene is, in general, stable and 
relatively good, individual nonlinearities and irregularities may be present. 
Accurate merging of satellite coordinates and latitude and longitude coordinates 
requires registry based on established ground-control points. Regression procedures 
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are used to develop equations relating the Landsat scene coordinates to 
latitude/longitudc coordinates at each control point. A detailed discussion of 
this procedure by Fisher and others is published in "Photograsunetric Eri^ . sring 
and Remote Sensing”, April 1979* Control points usually selected are conspicuous 
land features, such as penninsulas and islands. Residual errors of less than 
tvo pixels (picture element) in predicted versus actual row-column coordinates 
were determined. 

The basin boundary latitude and longitude coordinates stored in the cartographic 
base are converted, by the geometrically controlled regression equations, to 
Landsat row-column coordinates. The selected drainage basin is outlined, and 
the total amount of each land-use type within that individual basin is tabulated 
digitally and printed out (fig. 3). 

Landsat coordinates were determined for an area on the west side of Madison, 
Wis., to verify the? accuracy of the conversion technique. Land use was also 
classified to verify the classification technique. The predicted basin boundary 
agreed, to actual location with an accm-acy error of +_ one pixel. The land-use 
classifications agreed with visual spot checks. 


LOW-FLOW REGRESSION ANALYSIS 

The land-use data obtained from the Landsat imagery were used in a step- 
forward regression analyses, as outlined by Thomas and Benson ( 1970 ), to develop 
relations for low flow and floodflows. The equations with the lowest standard 
error of estimate with all variables significant at the 99 percent or the 95 percent 
confidence level were selected as the best equations for prediction. The resultant 
equations were compared to those developed using conventional land-use data. 

Low-flow relations were developed for the lower Wisconsin River basin 
(fig. U) based on data for 27 subbasins. The standard error of estimate (S£) 
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C5436250 


05436250 
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Part of a 1 :24,000 U.S.G.S. topographic quadrangle 
used to digitize the Allen Creek near Albany basin. 


Digitizcvi output at 1:250.000 
of Allen Creek near Albany basin. 



Tabular total for Allen Cr near Albany 
No. Theme No. of Pixels 


Unclassified 

Forest 


Wetlands 
Mixed veg. 
Bare toil 
Crop land 


Land use classification of Allen Creek 
near Albany. 


Figure 3. Steps for extracting land-use areas for a basin. 







fcr the relation g (leaet annual 7-day low flow for 2-year recurrence intervals) 
developed using Landsat iiaagcry was 35 percent (table l) as compared to an SE of 
1<1; percent fcr the relation developed using conventional land-use data. . .e SE 
for the ( 10-year recurrence intervals) relation was reduced from 52 io 

l»l percent (table l) when Landsat data were used. A reduced standard error of 
estimate indicates less variance about the relation, therefore the predictive 
capability of equation is improved. 

Gebert (1977) showed that the Q.^ g and relations were improved by 

using an additional basin characteristic, the base-flow index (Gebert, 1977. 
p. 15). The improved relations reported by Gebert are 'quations 3 and 5 in 
table 1. Through the use of iund-use data obtained from Landsat Imagery and 
including the base-flow index, the DE for the g relation decreased from 23 to 
19 percent, and the SE for the relation decreased from 32 to 26 percent 

fclc / (table l). 

Land-use data based on Landsat imagery for b3 subbasins (fig. b ) were used 
to develop Q,, „ and 0 relations for the Pecatonica-Sugar River basin as 

well. The SE's for these new relations showed no reduction. This occurred 
because only area and the base-flow index were significant i.-. the relations. 

Tlie land use of this basin apparently does not affect low flow as it does in the 
lower Wisconsin River basin. 

FLOOD-l’RHQUEI.'CY REGRESSION ANALYSIS 

Flood-frequency relations were also developed for the sections of the lower 
Wisconsin River and the Pecatonica-Sugar River basins in the "Driftless Area". 

The relations were developed through the use of Landsat land-use data, and these 
were compared to similar relations developed from conventional sources. Data 
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Table 2.— Comparison of flood-frequency regression resxilts through the use 
of Landaat data In the "Driftless Area" of Wisconsin 
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for lU subbasins were used to develop relations representing flood peaks vith 
recurrence intervals of 10, 50, and 100 years P^q» and P^qq)- The and 

P^Q relations developed through the use of Landsat imagery had SE's of aid 
17 percent, respectively, compared vith SE's of 29 and 31 percent, respectively, 
for the relations developed using conventional land-use data (table 2). The SE 
for the P^QQ relation was reduced from 3*+ to 17 percent when Landsat data were 
used. These represent significant improvements in the predictive relations. 

SUMMARY 

The Q.^ 2 0^ 3^0 Tow-flow relations in the lower Wisconsin River basin 

were improved 17 to 20 percent when Landsat land-use data were used in developing 
the relations. Land use did not prove to be a significant variable in the 
Pecatonica-Sugar River basin for low- flow estimates. Estimates for the flood- 
frequency relations with recurrence intervals of 10, 50 and 100 years were 
improved 1*5 to 50 percent. 

Lhndsat data can be used to significantly improve relations for predicting 
strearcflow characteristics. These improvements will result in a more accurate 
determination of water resources and will aid in management decisions. 
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ILLUSTRATIONS 


Figure 1. Map shoving location of the Pecatonica-Sugar and lover 
Wisconsin River basins in Wisconsin 
Figure 2. Annotated color-conposite Landsat-1 image of southvestern 
V7isconsin (5386-15325) 

Figure 3- Steps for extracting land-use areas for a basin 
Figure U. Map showing location of basins used for estimating low-flow 
characteristics 



